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1 Introduction 



Since their discovery in the mid-1970s, heavy-quarkonium states have played an im- 
portant role in helping us to understand QCD. Because heavy-quarkonium states 
are nonrelativistic, they allow the application of theoretical tools that simplify and 
constrain the analyses of nonperturbative effects. Hence, heavy-quarkonium states 
provide a unique laboratory in which to explore the interplay between perturbative 
and nonperturbative effects in QCD. 

In this talk I begin by describing the theoretical framework for calculations of 
inclusive quarkonium production processes, with emphasis on the NRQCD factor- 
ization approach [1]. I then describe the current status of the comparison between 
theoretical predictions and experimental measurements.* 

2 Theoretical Framework 

A number of theoretical approaches have been proposed for the calculation of heavy- 
quarkonium production processes. These include the NRQCD factorization approach 
[1], the fragmentation approach [4, 5], the color-singlet model (CSM) [6, 7, 8, 9, 10], 
the color-evaporation model (CEM) [11, 12, 13, 14], and the /c^-factorization approach 
[15, 16, 17, 18]. Only NRQCD factorization and the fragmentation approach are 
believed to be methods that could be derived from QCD. NRQCD factorization is 
the default model for most current studies of quarkonium production. 

2.1 NRQCD Factorization of the Inclusive Production Cross 
Section 

Nonrelativistic QCD (NRQCD) is an effective field theory that describes the behavior 
of bound states of a heavy-quark (Q) and a heavy-antiquark(Q) when the velocity v of 
the Q or Q in the QQ rest frame is nonrelativistic (d < 1). Some years ago, Bodwin, 
Braaten and Lepage [1] conjectured that the inclusive cross section for producing 
a quarkonium at large momentum transfer (p? or p*) can be written as a sum of 
products of "short-distance" coefficients times NRQCD matrix elements: 

*(!0 = 5>»(A)<0|O?(A)|0). (1) 

n 

The "short-distance" coefficients F n (A) are essentially the pro cess- dependent partonic 
cross sections to make a QQ pair, convolved with the parton distributions of the in- 
coming hadrons. The NRQCD long-distance matrix elements (LDMEs) (0|£^f(A)|0) 

*See Rcfs. [2, 3] for some further details regarding these issues. 
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are the probability for a QQ pair to evolve into a heavy quarkonium. A is the fac- 
torization scale, which is the cutoff of the effective field theory. The NRQCD matrix 
elements are vacuum expectation values of four-fermion operators in NRQCD, but 
with a projection onto an intermediate state of the quarkonium H plus anything: 

Of (A) = (0| X W \H + X)(H + X\ V<x|0>. (2) 
^ x ' 

Here, ^ and x are two-component (Pauli) fields that create a heavy quark and a 
heavy antiquark, respectively, and n n and n' n are direct products of Pauli and color 
matrices. ^ 

The short- distance coefficients have expansions in powers of a s . The LDMEs are 
nonperturbative, but they are conjectured to be universal, i.e., process independent. 
Only the color-singlet production LDMEs are simply related to the decay LDMEs. 
The LDMEs have a known scaling with v, where v 2 ~ 0.23 for the J/ip and v 2 ps 0.1 
for the T(1S). Hence, the NRQCD factorization formula is a double expansion in 
powers of a s and v. 

The current phenomenology of J/ip, ip(2S), and T(nS) production uses LDMEs 
through relative order w 4 : 



(O h ( 3 S?)) (O(v )), 

(o H esf ] )) (o(v*)), 

(O h ( 3 S?)) (0(v% 

(0 H ( 3 Pf ] )) (0(v% 



Here, the superscripts H stand for any spin-triplet, S'-wave quarkonium state, includ- 
ing the J/ip, the ijj(2S), and the T(nS) states, all of which have different NRQCD ma- 
trix elements. The quantities in parentheses are the spin, orbital-angular-momentum, 
and color (singlet or octet) quantum numbers of the QQ pair that evolves into the 
quarkonium state H. 

A key feature of NRQCD factorization is that quarkonium production can occur 
through color-octet, as well as color- singlet, QQ states. The color-singlet production 
LDMEs (O^S! 11 )} can be determined from quarkonium electromagnetic decay rates 
(up to corrections of order v 4 ). However, the color-octet LDMEs must be determined 
through comparisons of theoretical predictions with measurements. Once they have 
been fixed through comparisons between theory and experiment in one or more pro- 
cesses, the LDMEs can be used, by virtue of the property of universality, to make 
predictions for other processes. 



'It was pointed out by Nayak, Qiu, and Stcrman that gauge invariance requires that the defini- 
tions of the NRQCD LDMEs include Wilson lines that run from the quark and antiquark fields to 
infinity. For simplicity, I have omitted these Wilson lines here. 
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One obtains the CSM by retaining, for a given process, only the contribution that 
is associated with the color-singlet LDME of the lowest nontrivial order in v. The 
CSM is theoretically inconsistent in that it leads to uncanceled infrared divergences 
in calculations of production and decay processes for P-wave quarkonium states. 

2.1.1 Status of a Proof of NRQCD Factorization 

A proof of NRQCD factorization is complicated because gluons can dress the quarko- 
nium basic production process in ways that apparently violate factorization in indi- 
vidual Feynman diagrams. A proof of factorization would involve a demonstration 
that diagrams in each order in a s can be re-organized so that (1) all soft singularities 
cancel or can be absorbed into NRQCD LDMEs and (2) all collinear singularities 
and spectator interactions can be absorbed into parton distributions. Nayak Qiu 
and Sterman have demonstrated factorization for the leading power in 1 /p? through 
next-to-next-to- leading order (NNLO) in a s . However, it is not known if the proof 
generalizes to all orders in a s . Note that an all-orders proof is essential because 
potential violations of factorization involve soft gluons, for which a s is not a good 
expansion parameter. It seems likely that, if a proof of NRQCD factorization can be 
established, it would hold only for values of pt that are greater than the heavy-quark 
mass rriQ. 

2.2 The Problem of Large Higher-Order Corrections 

2.2.1 Higher-Order Corrections in J/ip and T(1S) Production 

When one calculates quarkonium production cross sections differential in p? in the 
NRQCD factorization approach, large corrections appear in higher orders in pertur- 
bation theory. These large corrections were first noticed in calculations of quarkonium 
production at hadron-hadron colliders through the color-singlet channel [19]. An ex- 
ample of this phenomenon is shown in Fig. 1. The curve labeled "NNLO*" is an 
estimate that is based on only the real-emission contributions at NNLO in a s . Note 
that the color-singlet contributions alone do not explain the data. 

An example of large higher-order corrections in T(IS') production in the color- 
singlet channel is shown in Fig. 2. The next-to-leading order (NLO) results that are 
shown in this figure were confirmed by Gong and Wang [22]. In this case, the data 
could be explained by NNLO* color-singlet production alone. However, given the 
large error bars, there is still room for a substantial of color-octet contribution. 

A very large correction in next-to-leading order in a s also arises in J / ifj production 
through the color-octet 3 Pj channel, where a negative correction has been found 
[23, 24]. However, NLO corrections to the S-w&ve color-octet channels are small [25]. 
The corrections for production at the Tevatron are approximately factors of 1.235 for 
the l So color-octet channel and 1.139 for the 3 Si color-octet channel. 
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Figure 1: Calculations of the J/ip production cross section in the color-singlet channel 
compared to CDF Run II data [20]. LO: leading order in a s . NLO: next-to-leading- 
order in a s . NNLO*: next-to-next-to-leading order in a s , real corrections only. (Fig- 
ure from P. Artoisenet, based on work by P. Artoisenet, J. Campbell, J. P. Lansberg, 
F. Maltoni, and F. Tramontano [19].) 




Figure 2: Calculations of the T(1S) production cross section in the color-singlet 
channel compared to CDF Run II data [21]. LO: leading order in a s . NLO: next- 
to-leading-order in a s . NNLO*: next-to-next-to-leading order in a s , real corrections 
only. (Figure from P. Artoisenet, based on work by P. Artoisenet, J. Campbell, 
J. P. Lansberg, F. Maltoni, and F. Tramontano Ref. [19].) 

These results raise several questions: (1) Does the perturbation series converge? 
(2) How does one understand the different sizes of the higher-order corrections for 
different channels? 
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2.2.2 Explanation of the Large Corrections 

It has been suggested by Campbell, Maltoni, and Tramontano [26] that the large 
higher-order corrections to quarkonium production can be understood as follows: At 
high px, higher powers of a s can be offset by a less rapid fall-off with p?. For example, 
for production of the J ftp through the S-w&ve, color-singlet channel, the contribution 
at leading order (LO) in a s goes as a 3 rn^/p\,. In contrast, the NLO contribution goes 
as a^m^/p^ and is enhanced at high p? relative to the LO contribution. * The NNLO 
contribution to S- wave, color-singlet production of the J/ip is further enhanced at 
high px, since it goes as a^/p^, owing to contributions in which a gluon fragments 
into the J/ip. The leading power behavior of the cross section at high px is 
[4, 5]. Hence, one expects no further large kinematic enhancements beyond NNLO. 

The NNLO calculation of color- singlet, S"-wave production of S- wave charmo- 
nium states is, so far, incomplete, as only the real-emission (NNLO*) correction has 
been computed. It has been noticed by Ma, Wang, Chao [27] that the color-singlet 
NNLO* correction seems be dominated by contributions that are proportional to 
log 2 \p\ I ^T CUt ) , where p\ mt > trq is an artificial cutoff that is introduced to make 
the calculation finite. These contributions will cancel when virtual corrections are 
included in the complete NNLO calculation. Hence, the complete NNLO contribu- 
tion is likely to be significantly smaller than the NNLO* contribution and cannot, by 
itself, explain the observed production cross sections of S"-wave charmonium states. 

Kinematic enhancements also account for the large corrections in the color-octet 
channels. The color-octet 3 S\ channel receives a small correction in NLO because 
it already has 1/p? behavior in LO, owing to a gluon-fragmentation contribution. 
The color-octet 3 Pj channel receives a large correction in NLO because it first shows 
1/pr behavior in NLO, again owing to gluon fragmentation. The color-octet 1 So 
channel also shows first shows \jp% behavior in NLO, and that behavior arises from 
gluon fragmentation. However, the NLO correction to the color-octet 1 So channel is 
numerically small at moderate pr because the fragmentation process has little support 
near z = 1} 

Note that, at NLO, all of the color-octet channels that appear in 5-wave quarko- 
nium production already contain contributions that go as l/p?- This suggests that 
no further large enhancements will occur in still higher orders in a s in these channels. 
At leading-order (LO) in a s , only the 3 Si color-octet channel has a behavior. 
This fact has been used to argue that the 3 Si color-octet channel is dominant at large 

^The NLO correction also contains a contribution, which arises from charm-quark fragmentation 
into a J/ip, that goes as ct^/pi^. This charm-quark fragmentation contribution happens to be small 
numerically. 

§The absence of peaking near z = 1 may reflect the fact that the color-octet 1 So fragmentation 
process contains no soft divergences in full QCD [28]. Such soft divergences are ultimately absorbed 
into the NRQCD LDMEs, rendering the physical cross section finite, but the remnants of the soft 
divergences can produce a peaking near 2 = 1. 
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Pt- However, at NLO, all of the color-octet channels have a l/p? behavior, and so the 
argument for the dominance of the 3 Si color-octet channel at large px is not valid. 

2.3 The Fragmentation Approach 

Kang, Qiu, and Sterman [4, 5] have suggested that one reorganize the perturbation 
expansions for inclusive quarkonium production cross sections according to the Pt 
dependence of the various contributions. Specifically, they find that the leading be- 
havior in 1/pt in the production cross section (1/pj-) comes from contributions in 
which single-parton production cross sections are convolved with the fragmentation 
functions for a single parton into a quarkonium: 

d&A+B^i+x <S> A-Kff- (4) 

They also find that the first subleading behavior in 1/pr in the production cross sec- 
tion {itiq/pt) comes from QQ production cross sections convolved with fragmentation 
functions for a QQ pair into a quarkonium: 

d&A+B^QQ+X ® D Q q^ h . (5) 

Nayak, Qiu, and Sterman have provided convincing arguments that these results hold 
to all orders in perturbation theory, up to corrections of order ttlq/pj.. 

The fragmentation approach of Kang, Qiu, and Sterman is well suited to the 
analysis of phenomenon of large higher-order corrections because it allows one to 
identify the specific sources of the kinematic enhancements at high p^, which arise 
from the leading and first subleading powers in 1/pt in the quarkonium production 
rates. Hence, the fragmentation approach should allow one to reduce the uncertainties 
in the theoretical predictions by focusing on those processes that produce the large 
corrections, calculating them to yet higher orders in a s , and resumming the associated 
logarithms of p^/rriQ, which are large at high p?. 

It is important to check that the fragmentation formalism really does account for 
all of the large corrections. It has been confirmed that the fragmentation contribution 
reproduces most of the large correction in the color-singlet channel at NLO [5]. There 
is also an estimate of the fragmentation contribution in the 3 Pj color-octet channel 
that indicates that it gives a significant part of the NLO correction in that channel 
[28]. Several other tests of the fragmentation approach are in progress. 

The values of the nonperturbative fragmentation functions that appear in this 
formalism are, so far, unknown. However, if NRQCD factorization holds, then the 
fragmentation functions of Kang, Qiu, and Sterman can be written as a sum of 
NRQCD LDMEs times perturbatively calculable short- distance coefficients. In that 
case, the fragmentation approach is equivalent to a reorganization of the NRQCD 
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factorization formula (1) according to the behaviors and tuq/p^,, with the stip- 
ulation that there are corrections of order frig/p? that are not accounted for in the 
fragmentation approach. 

3 Comparisons of NRQCD Factorization with Ex- 
periment 

3.1 Overview 

NLO corrections to charmonium production cross sections and polarizations have been 
computed for many production processes: J/ip and ip(2S) production cross sections 
and polarization at the Tevatron and the LHC; J /if) and if>(2S) production cross 
sections at RHIC; J /if) photoproduction cross sections and polarization at HERA; 
the J/ip + r] c production cross section, the J /if) + cc production cross section, and the 
J/if) + X(non-cc) production cross section in e + e~ annihilation at the B factories. An 
NLO calculation of the xj production cross sections for the individual J states also 
exists [29]. 

Generally, data and the NLO predictions of NRQCD factorization for quarkonium 
production agree, within errors.^ There are three significant exceptions: (1) quarko- 
nium polarization in hadron-hadron collisions , (2) the cross section for the production 
of J /if) + X(non-cc) in e + e~ annihilation at Belle, and (3) the cross section for J /if) 
production in 77 scattering at LEP II. I will discuss each of these exceptions below, 
following a discussion of the extraction of the NRQCD LDMEs from hadroproduction 
and photoproduction data. 

3.2 Extraction of NRQCD LDMEs at NLO 

I now describe an important issue that arises in extracting NRQCD matrix elements 
from fits of NLO NRQCD predictions to J /if) production data. 

Recently, two groups, Ma, Wang, and Chao [23, 27] and Butenschon and Kniehl 
[24], have carried out the first complete calculations of J /if) hadroproduction at NLO. 
These calculations are complete in the sense that they include all of the color-singlet 
and color-octet channels that contribute through order v 4 . [See Eq. (3).] The results 
of the two groups for the NRQCD short- distance coefficients agree. However, the 
fitted NRQCD matrix elements are very different. 

Using the CDF Run II data, Ma, Wang, and Chao [23, 27] could fit only two linear 

^See the global fit of NRQCD matrix elements of Butenschon and Kniehl [30, 31] and Rcf. [3] for 
some further details. 
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combinations of matrix elements unambiguously: 

M , ro = (O*^)) + (r /m 2 c )(O^CP?)) = (7.4 ± 1.9) x 1(T 2 GeV 3 , 

M 1>ri = (0^( 3 5f)) + {n/m 2 c )(O^^Pf ] )) = (0.05 ± 0.03) x 10~ 2 GeV 3 , (6) 

where ro = 3.9 and r\ = —0.56 were chosen on the basis of approximate relations 
between the short-distance coefficients. 

Butenschon and Kniehl [30] used their NLO calculations for ep, 77, and e + e~ 
production to fit all three color-octet LDMEs, using data from the Tevatron, LHC, 
RHIC, HERA, LEPII and KERB: 

(o^4 8] )> 

which implies that 

M , ro = (2.17 ± 0.56) x 10" 2 GeV 3 , 

M ljn = (0.62 ± 0.08) x 10" 2 GeV 3 . (8) 

There are many small differences in the fitting procedures: (1) Ma, Wang, and 
Chao included feeddown from the i/j{2S) and the Xcj states in the theoretical predic- 
tion, but Butenschon and Rniehl did not; (2) Ma, Wang, and Chao used a 2-parameter 
constrained fit, while Butenschon and Rniehl used a 3-parameter fit; (3) Ma, Wang, 
and Chao applied a cut px > 7 GeV to the CDF data, while Butenschon and Rniehl 
applied a cut pt > 3 GeV to the CDF data. However, the most important difference 
between these fits is the use of the 1ow-j?t HI data by Butenschon and Rniehl. Indeed, 
the fit of Butenschon and Rniehl to the Tevatron and HERA data alone [24] gives a 
very similar result to that of the global fit: 

M 0iro = (2.5 ± 0.08) x 10^ 2 GeV 3 , 

M 1>ri = (0.59 ± 0.02) x 10~ 2 GeV 3 . (9) 

Note that, for the HERA data, px lies in the range 1 GeV < px < 3 GeV. It is not 
clear that NRQCD factorization, if proven, would hold at such low values of px- 

Both fits describe the data within errors, but, as can be seen from Fig. 3, the shape 
of the Ma, Wang, and Chao fit agrees with the CDF data better than the shape of 
the Butenschon and Rniehl fit. The fit of Butenschon and Rniehl to the CDF data 
clearly overshoots the data at large px, while the fit of Ma, Wang, and Chao does not. 
Resummation of large logarithms of px/ml at large px would lower the theoretical 
predictions because the fragmentation process softens the J/tp momentum. This 



= (4.76 ± 0.06) x 10" 2 GeV 3 , 

= (0.265 ± 0.014) x 10" 2 GeV 3 , 

= (-0.716 ± 0.089) x 10~ 2 GeV 3 , (7) 
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Figure 3: Comparisons of NLO fits to the CDF data [32]. Left: The fit of Ma, Wang, 
and Chao [23, 27]. (Figure from Ref. [23].) Right: The global fit of Butenschon and 
Kniehl [30]. In this, and subsequent, figures "CS" is the color-singlet contribution 
and "CO" is the color-octet contribution. (Figure from Ref. [30].) 

resummation correction would improve the agreement of the fit of Butenschon and 
Kniehl with the CDF data and worsen the agreement of the fit of Ma, Wang, and 
Chao with the CDF data. Note that, in this figure (and in subsequent figures from 
the works of Butenschon and Kniehl), it is clear that the color-singlet contribution 
alone does not describe the data. 

As can also be seen from Fig. 3, the prediction of Ma, Wang, and Chao deviates 
from the data at low px- Presumably, the prediction of Butenschon and Kniehl would 
show a similar deviation if it were extended to sufficiently low Pt- This deviation may 
be an indication that factorization is breaking down, and/or that large logarithms of 
nrl/p T must be be resummed in the low-pr region. As can be seen in Fig. 4, there is 
also a slight discrepancy in shape in the fit of Butenschon and Kniehl to the HI data, 
which might be attributable to a breakdown of factorization and/or a need to resum 
logarithms of m 2 c /p T in the low-pr region. 

The fit of Ma, Wang, and Chao to the Tevatron data and the fit of Butenschon 
and Kniehl [24] to the Tevatron and HERA data have been used to make predictions 
for J/ijj cross sections at the LHC. Comparisons to the LHCb data to the prediction 
of Ma, Wang, and Chao and to the prediction of Butenschon and Kniehl are shown 
in Fig. 5. 

As can be seen, the prediction of Ma, Wang, and Chao yields a good fit to the 
data. The prediction of Butenschon and Kniehl also fits the data within errors, but 
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Figure 4: The global fit of Butenschon and Kniehl [30] compared to the HI data 
[33, 34]. (Figure from Ref. [30].) 
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Figure 5: Comparisons of NLO predictions with the LHCb data [35]. Left: The 
prediction of Ma, Wang, and Chao. (Figure from Ref. [35].) Right: The prediction 
of Butenschon and Kniehl from a fit to the Tevatron and HERA data [24]. (Figure 
from Ref. [35].) 



has the wrong shape: It overshoots the data at high pr- This shape discrepancy at 
high pr is even more apparent in Fig. 6, which shows a prediction of Butenschon 
and Kniehl, which is based of their global fit, in comparison to Atlas data. Again, a 
resummation of large logarithms of Px/ m c & t high pt may improve the shape of the 
Butenschon and Kniehl predictions, while worsening the shape of the Ma, Wang, and 
Chao predictions. 
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Figure 6: Comparison of a prediction based on the global fit of Butenschon and Kniehl 
[30] with Atlas data [36]. (Figure from Ref. [30].) 

Comparisons of PHENIX data with the prediction of Ma, Wang, and Chao and 
the global fit of Butenschon and Kniehl are shown in Fig. 7. The prediction of Ma, 




Figure 7: Comparisons of NLO calculations with PHENIX data [37]. Left: The 
prediction of Ma, Wang, and Chao [27]. (Figure from Ref. [27].) Right: The global 
fit of Butenschon and Kniehl [24]. (Figure from Ref. [24].) 

Wang, and Chao agrees well with the data. This prediction takes into account the 
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effects of feeddown from the Xcj states and the i/j{2S) in both the fit to the CDF data 
that is used to extract the LDMEs and in the prediction for the PHENIX data. The 
global fit of Butenschon and Kniehl undershoots the PHENIX data slightly, but still 
agrees within errors. Feeddown effects are not taken into account in the global fit of 
Butenschon and Kniehl, but they are expected to be small [24]. 

After this talk was presented, a new NLO calculation of J / hadropro duct ion by 
Gong, Wan, Wang, and Zhang [38] appeared. This calculation includes contributions 
at NLO from feeddown from the Xcj and ip{2S) states and extracts all three color- 
octet NRQCD LDMEs from a combined fit to the CDF Run II [20] and LHCb [35] 
measurements of the cross section, differential in pr- The fit makes use only of data 
with pt > 7 GeV. Gong et al. were able to determine three linear combinations of 
color-octet LDMEs with good accuracy. Their LDMEs yield 

M , ro = (6.0 ± 1.0) x 10~ 2 GeV 3 , 

Mi, n = (0.07 ± 0.01) x 10" 2 GeV 3 , (10) 
in good agreement with the fit of Ma, Wang, and Chao [Eq. (6)]. 

3.3 Polarization 

Polarization of a 3 S\ quarkonium state is analyzed experimentally in terms of the 
angular distribution the lepton pairs into which it decays. The angular distribution 
is usually parametrized as 

W(6,4>) oc 1 + A e cos 2 fl + A^sin 2 #cos(20) + A^sin(2#)cos0, (11) 

where 9 and <fi are the polar and azimuthal angles, respectively, of the positively 
charged lepton with respect to some coordinate system (often called the "frame") 
and Xg, A^, and A^ are the measured polarization parameters. The parameter A^, 
which is also called a in the literature, is a measure of the amount of transverse or 
longitudinal polarization of the quarkonium: A^ = 1 corresponds to completely trans- 
verse polarization, and A^ = —1 corresponds to completely longitudinal polarization. 

Quarkonium polarization at large px has long been thought to provide a definitive 
test of the color-octet production mechanism. However, this is an idea that is based 
on arguments at LO in a s , and, as we will see, it may not be correct at NLO. 

3.4 Polarization at LO in ot s 

At LO in a s , hadroproduction of a 3 Si quarkonium at large px is dominated by the 
process of gluon fragmentation into a color-octet 3 Si QQ pair. At large pr in the 
hadronic center-of-momentum (CM) frame, the gluon is nearly on mass shell, and, 
so, is transversely polarized. At LO in a s , that transverse polarization is transferred 
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completely to the QQ pair in the perturbative fragmentation process. As the QQ 
pair evolves nonperturbatively into a 3 Si quarkonium through the emission of gluons 
that are soft in the QQ rest frame, the transverse polarization is preserved at LO 
in v because of the heavy-quark spin symmetry [39]." The transverse polarization is 
defined in the hadronic CM frame relative to the direction of the fragmenting gluon's 
momentum, and is invariant under a boost in that direction. At LO in a s and v, that 
boost vector coincides with the boost vector from the quarkonium rest frame to the 
hadronic CM frame, which defines the polar axis of the "helicity frame." Hence, at 
LO in a s and v , one expects Xg in the helicity frame to approach 1 as pt increases. 

The LO prediction of transverse polarization at large px has not been confirmed 
experimentally. As can be seen from the left-hand plot in Fig. 8, the LO prediction 




20 p T (GeV) p T (GeV/c) 



Figure 8: Comparisons of the LO NRQCD prediction for polarization of the J/tp [41] 
with CDF data. Left: The CDF Run I data [42]. (Figure provided by Jungil Lee.) 
Right: The CDF Run II data [20]. (Figure from Ref. [20].) 

agrees, within errors, with the CDF Run I measurement, except for the two highest 
Pt points. However, as can be seen from right-hand plot in Fig. 8, the LO prediction 
disagrees completely with the CDF Run II measurement. Feeddown contributions 
from the Xc states (about 30%) and the ip(2S) state (about 10%) are included in both 
the data and the theoretical prediction. The CDF Run I and Run II measurements 
are also in disagreement with each other. The reason for this discrepancy has not 
been found. 

The LO NRQCD factorization prediction for polarization of the T(IS') is compared 
with the CDF and DO Run II data in Fig. 9. Both the CDF and DO results are 
incompatible with the LO prediction, but in different regions of px- The CDF and 

"There are spin- flip interactions, which are suppressed as v 3 . This suppression has been verified 
in a lattice calculation of decay matrix elements [40] . 
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Figure 9: Comparison of the LO NRQCD prediction for polarization of the T [41] 
with the CDF Run II data [43] and the DO Run II data [44]. (Figure provided by Hee 
Sok Chung, using Refs. [43, 44].) 

DO results are also incompatible with each other. No reason for the incompatibility 
of the experimental results has been found. 

3.4.1 First Complete NLO Calculations of Polarization 

The first NLO calculation of J/ip polarization in photoproduction was carried out 
by Butenschon and Kniehl [45]. This calculation is complete in the sense that it 
includes all of the production channels that appear in NRQCD through relative order 
v A [Eq. (3)]. It was preceded by a calculation at NLO for the 3 Si color-singlet channel 
for J/ip photoproduction. [46]. 

The first complete NLO calculations of J/ip polarization in hadropro duct ion were 
carried out by Butenschon and Kniehl [47] and by Chao, Ma, Shao, Wang and Zhang 
[48]. These calculations also include all of the production channels that appear in 
NRQCD through relative order v A . They were preceded by partial calculations at NLO 
for the 3 Si color-singlet channel for J/ip hadroproduction [22] and Y(1S) hadropro- 
duction [19] and by calculations at NLO for the 3 Si and 1 5'o color-octet channels for 
Jjijj hadroproduction [25] and Y(1S) hadroproduction [49]. 

At NLO, the theoretical predictions for polarization change dramatically. Polar- 
ization in the color-singlet channel changes from transverse to longitudinal [19, 22]. 
However, this channel makes a very small contribution to the hadroproduction cross 
section. A much more important issue is that, at NLO at high p?, there are large 
corrections to production in the 3 Pj color-octet channel that are mostly transversely 
polarized. At NLO, the 3 Si color-octet channel is still transversely polarized at high 
Pt- However, owing to the large NLO corrections in the 3 Pj color-octet channel, it 
is no longer the dominant channel for hadroproduction at high p?, and it is not the 
only source of transverse polarization. 
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Fig. 10 shows the prediction of Butenschon and Kniehl for J/ip polarization in 
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Figure 10: Comparisons of the NLO NRQCD prediction of Butenschon and Kniehl 
[47] for the polarization of the J/ip with hadroproduction data. Left: the CDF Run I 
[42] and Run II [20] data. (Figure from Ref. [47].) Right: Alice data [50]. (Figure 
from Ref. [47].) 

hadroproduction in comparison with the CDF Run I and Run II data and Alice data. 
The prediction is based on the LDMEs from the global fit of Butenschon and Kniehl 
[30]. It is in considerable disagreement with the data. On the other hand, as can be 
seen from Fig. 10, the NLO J/ip polarization prediction of Butenschon and Kniehl is 
in rough agreement with the Alice data. It should be noted, though, that the Alice 
data include J/ip production from 5-meson decays and feeddown from the Xcj and 
i(>(2S) states. 

The LDMEs from the global fit of Butenschon and Kniehl can also be used to 
predict the polarization of the J /if) in photoproduction at HERA. That prediction is 
compared in Fig. 11 with HI data. The data are roughly compatible with the theory 
at large p?, but the error bars are large. 

In contrast, Chao, Ma, Shao, Wang, and Zhang [48] determined the three color- 
octet NRQCD LDMEs by using data from the CDF Run II measurements of the J ftp 
cross section, differential in px, and the CDF Run I and Run II measurements of 
the J /ip polarization, differential in p?. The resulting LDMEs are barely compatible 
with the LDMEs from the fit of Ma, Wang, and Chao [23, 27] to the production 
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Figure 11: Comparison of the NLO NRQCD prediction of Butenschon and Kniehl 
[47] for the polarization of the J/ip with HI data [34, 51]. (Figure from Ref. [47].) 

cross section alone: The values of M 0jro and M l ri that Chao et al. obtain differ from 
those of Ma, Wang, and Chao by about 2a. This tension between the two sets of 
LDMEs suggests that the LDMEs from the fit of Ma, Wang, and Chao might lead to 
a polarization prediction that deviates from the CDF Run II data, although perhaps 
with large uncertainties. 

The LDMEs of Chao et al. lead to the prediction for J/ip polarization at CDF 
that is shown in Fig. 12. With this choice of LDMEs, the contributions of the color- 
octet 3 Si and 3 Pj channels to the transverse polarization largely cancel, leaving a 
net polarization parameter \q that is near zero, in agreement with the CDF Run II 
data. This choice of LDMEs still gives reasonable predictions for the J/ip production 
cross sections that have been measured by the Atlas and CMS Collaborations, as 
can be seen from Fig. 13. This reflects the fact that the J/ip LDMEs and, hence, 
the J/ip polarization in hadroproduction, are not well determined from fits to the 
production cross section alone. The Chao et al. LDMEs lead to a prediction of a 
slight longitudinal polarization at the LHC, as can be seen in Fig. 14. However, in 
Fig. 15, one can see that the Chao et al. LDMEs seem to be incompatible with the 
HERA photoproduction data, even at large px- 

After this talk was presented, Gong et al. also completed a new NLO calculation 
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Figure 12: Comparison of the NLO NRQCD prediction of Chao, Ma, Shao, Wang, 
and Zhang [48] for the polarization of the J/ip with the CDF Run I [42] and Run II 
[20] data. (Figure from Ref. [48].) 




Figure 13: Comparisons of LHC data with the NLO NRQCD prediction for the J/tp 
cross section of Chao, Ma, Shao, Wang, and Zhang [48]. The prediction uses NRQCD 
LDMEs that are extracted from fits to the CDF measurements of the J/tp cross section 
[32] and the J/ip polarization [20, 42]. Left: The Atlas data [36]. Right: The CMS 
data [52]. (Figure from Ref. [48].) 

of J/ip polarization in hadroproduction [38], including the effects at NLO of the 
feeddown contributions from the Xcj and ip(2S) states. Using the color-octet LDMEs 
from their fits to the CDF Run II and LHCb data, Gong et al. made predictions for 
the J /ip polarization. Their prediction for the polarization at CDF is intermediate 
between that of Butenschon and Kniehl [47] and that of Chao et al. [48], agreeing 
with the CDF Run I data [42], except for the highest-pr point, but disagreeing with 
the CDF Run II data [20]. Since the polarization predictions of Gong et al. rely only 
on data with p? > 7 GeV, they are less subject to questions about the validity of 
factorization than the predictions of Butenschon and Kniehl [47]. 
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Figure 14: NLO NRQCD prediction of Chao, Ma, Shao, Wang, and Zhang [48] for 
the polarization of the J/ip at the LHC. (Figure from Ref. [48].) 



10 





60 GeV < W < 240 GeV : 
0.3<z<0.9 
5 ^\ Q 2 < 1 GeV 2 I 
S VI - 314 GeV : 




"•<:;;.. i ^ 






• CS, LO 

- CS, NLO 

- CS+CO, LO 
- CS+CO, NLO 


5 \ \ 

I \\ ] 




H1 data: HERA1 
i 




10 

P? [GeV 2 ] 



Figure 15: Comparison of a prediction based on the LDMEs of Chao, Ma, Shao, Wang, 
and Zhang [48] with the HI data [33, 34]. (Figure provided by Mathias Butenschon.) 

3.5 e + e~ — >> J/i(j + X(non-cc) 

The total cross section for the process e + e~ — > J ftp + X(non-cc) has been measured 
by the Belle Collaboration [53] . Their result is 

a(e + e~ -> J/ip + X(non-cc)) = 0.43 ± 0.09 ± 0.09 pb. (12) 

NLO calculations using the NRQCD factorization approach have been carried out by 
Zhang, Ma, Wang, and Chao [54] and by Butenschon and Kniehl [30] and lead to the 
prediction 

a(e + e~ -> J/ijj + X(non-cc)) = 0.99±g;?| pb (/i = y/s/2), (13) 
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which is based on the NRQCD LDMEs from the global fit of Butenschon and Kniehl 
[30] and includes the effects of feeddown from the \cj and ip(2S) states [54]. The 
comparison of the theoretical calculations with the measured cross section favors the 
value of Mo )ro from the global fit of Butenschon and Kniehl, rather than the value 
from the fit of Ma, Wang, and Chao [23, 27]. 

There is tension between the Belle data and the theoretical prediction. However, 
it is worth noting that the most recent Belle measurements [53] imply that 

a(e + e~ -> J/ip + X) = a(e + e~ -> J/ip + cc + X) + a(e + e~ ->• J/ip + X(non-cc)) 

= 1.17±0.12ig;l|pb, (14) 

while the BaBar Collaboration obtained [55] 

a(e + e~ -> J /if) + X) = 2.52 ± 0.21 ± 0.21 pb. (15) 

This suggests the possibility that the Belle value of <r(e + e~ — > J /ip + X(non-cc)) is 
too small. Furthermore, most of the Belle data are at px < 3 GeV. Even if NRQCD 
factorization could be established at high px, it might not hold at such small values 
of pt- 

3.6 J 1^ Production in 77 Scattering at LEP II 

Butenschon and Kniehl have also made an NLO calculation of J/ip production in 
77 collisions at LEPII and incorporated it into their global fit [30]. This global fit is 
compared with DELPHI data in Fig. 16. The DELPHI data are slightly incompatible 
with the global fit. However, the error bars are large, especially at large px, and one 
should keep in mind that NRQCD factorization, even if it is established at high px, 
may not hold at such low values of px- 

4 Conclusions 

Studies of quarkonium production provide a unique laboratory in which to understand 
the interplay between the perturbative and nonperturbative aspects of QCD. They 
also provide an important testing ground for techniques for perturbative calculations 
of high-energy production cross sections. 

At present, the primary method for theoretical calculations of quarkonium pro- 
duction is the NRQCD factorization approach. NRQCD factorization provides a 
systematic framework for computing quarkonium production. However, there is no 
proof of NRQCD factorization for quarkonium production beyond two loops. Fur- 
thermore, in spite of considerable theoretical and experimental progress over the last 
decade, a definitive test of NRQCD factorization has not yet been achieved. 
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Figure 16: Comparison of the global fit of Butenschon and Kniehl [30] with the 
DELPHI LEP II data [56]. (Figure from Ref. [30].) 

Large corrections at NLO in a s to inclusive quarkonium production are now be- 
lieved to be understood in terms of kinematic enhancements. This understanding 
suggests that, for the color-octet production channels, which are believed to be much 
more important than the color-singlet production channels, there will be no large 
kinematic enhancements beyond NLO. If that is the case, then the existing NLO cal- 
culations of quarkonium production should at least be qualitatively correct. However, 
it seems likely that resummation of large logarithms will be needed in order to bring 
theory into good agreement with experiment at the large values of quarkonium px 
that will be probed at the LHC. 

The predictions of the color-singlet model fail to describe the data. In contrast, 
the predictions of NRQCD factorization are in agreement with most of the inclusive 
production data. There are several notable exceptions: (1) J/ip and T(2S) polariza- 
tion at the Tevatron, (2) the e + e~ — >■ J/ip + X(non — cc) total cross section, and (3) 
J/ip production in 77 scattering at LEP II. There are experimental discrepancies that 
cast some doubt on the polarization measurements and on the measurement of the 
e + e~ — >■ J/ip + X(non — cc) total cross section. The measurement of J/ip production 
in 77 collisions has large experimental uncertainties. Furthermore, even if NRQCD 
factorization is established at high p?, it might not hold at the low values of pt at 
which this measurement was made. 

The J /if) cross sections measured at the Tevatron and HERA, the J/ip polariza- 
tion measured by the CDF Collaboration in Run II, and NRQCD factorization seem 
to be mutually incompatible. Furthermore, the new polarization prediction of Gong 
et al. [38], which appeared after this talk was presented, suggests that the J/ip cross 
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sections measured by the CDF and LHCb collaborations, the J /if) polarization mea- 
sured by the CDF Collaboration in Run II, and NRQCD factorization are mutually 
incompatible. Some possible explanations for these incompatibilities are: (1) the 
CDF Run II J /if; polarization measurement is wrong; (2) NRQCD factorization fails 
at the low values of pt at which the HERA measurements were made; (3) additional 
corrections to the theoretical predictions are needed, such as those from higher orders 
in v; (4) NRQCD factorization is incorrect. 

Further advances in the theory of quarkonium production would help to sort out 
the various possibilities. These include a proof or a disproof of NRQCD factoriza- 
tion, calculations of corrections of higher order in a s and/or resummations of large 
logarithms of p\/m 2 c for the fragmentation processes that are dominant at large Pt, 
and calculations of the rates of additional production processes that could help to pin 
down the NRQCD LDMEs and to provide further tests of the production mechanisms. 

Additional experimental measurements are essential to the process of gaining an 
understanding of the quarkonium production mechanisms. Experimental efforts that 
would be of immediate value are (1) measurements J /if; direct production cross sec- 
tions at high pt] (2) measurements of all three J/if> polarization parameters in differ- 
ent frames for direct production at high p?; (3) measurements of XcJ cross sections 
and polarizations at high p?; (4) measurements of T direct production cross sections 
at high pt] (5), measurements of additional high-p^ quarkonium production cross 
sections, such as the J /if) + jet cross section at high pr- 

Since NRQCD factorization, if it is correct, likely holds only for px substantially 
greater than tuq, it is essential that all measurements be made at the highest possible 
values of p?. It is also very important to make measurements of direct, rather than 
prompt, cross sections. The feeddown contributions in prompt cross sections greatly 
complicate the theoretical analyses, which may not be under good control even in 
the simpler case of direct production. The measurements of Xcj cross sections and 
polarizations might be particularly interesting because the theory is more constrained: 
Only two NRQCD LDMEs enter at the leading non-trivial order in v. T cross sections 
and polarizations provide another important test of NRQCD factorization because the 
lower value of v 2 in the T(IS'), in comparison with the J / if>, means that the T(IS') and 
J / if; LDMEs have different relative sizes. It is clear that measurements of quarkonium 
production processes and polarizations that take advantage of the high px reach and 
high luminosity of the LHC will play a large role in advances in our understanding of 
quarkonium production mechanisms in the next few years. 
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